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The BBN predictions for the abundances of the light element isotopes is reviewed 
and compared with recent observational data. The single unknown parameter of 
standard BBN is the baryon-to-photon ratio, n, and can be determined by the con- 
cordance between theory and observation. Recent CMB anisotropy measurements 
also lead to a determination of n and these results are contrasted with those from 
BBN. In addition, the CMB data indicate that the Universe is spatially fiat. Thus 
it is clear that some form of non-baryonic dark matter or dark energy is necessary. 
Here I will also review the current expectations for cold dark matter from minimal 
supersymmetric models. The viability of detecting supersymmetric dark matter 
will also be discussed. 



1 Introduction 

The cornerstones of the Big Bang theory are the cosmic microwave background 
radiation (CMB) and big bang nucleosynthesis (BBN) (one could make an 
argument to include inflation as well). The existence of the former and the 
success of the latter point unequivocally to a hot and dense origin to the 
Universe. Indeed, these two aspects of the theory are intimately linked as 
early work by Alpher and Hermana on BBN led to their prediction of the 
CMB with a temperature of order 10K. 

These lectures will focus on recent developments in BBN theory and the 
related observations which test for concordance. Recent measurements of the 
CMB power spectrum have allowed an unprecedented level of accuracy in the 
determination of cosmological parameters including the baryon density, which 
is the key parameter for BBN. The concordance between BBN and the recent 
CMB measurements will be addressed. It is also becoming clear that while 
the total density of the Universe is near critical, i.e., we live in a spatially 
fiat Universe, the baryon density is only a small fraction of the total energy 
density. Something is missing. While much of missing energy appears to be 
best fit by a smooth component such as a cosmological constant, a sizable 

♦SUMMARY OF INVITED LECTURES AT THE FIRST NCTS WORKSHOP ON 
ASTRO PARTICLE PHYSICS, DECEMBER 6 - 9, 2001, RENTING, TAIWAN. 
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fraction must lie in the form of non-baryonic dark matter. The second half 
of these lectures will focus on the specific possibility of supersymmetric dark 
matter. 

2 Big Bang Nucleosynthesis 

The standard modelH of big bang nucleosynthesis (BBN) is based on the rel- 
atively simple idea of including an extended nuclear network into a homoge- 
neous and isotropic cosmology. Apart from the input nuclear cross sections, 
the theory contains only a single parameter, namely the baryon-to-photon 
ratio, 77. Other factors, such as the uncertainties in reaction rates, and the 
neutron mean-life can be treated by standard statistical and Monte Carlo 
techniques! The theory then allows one to make predictions (with well- 
defined uncertainties) of the abundances of the light elements, D, 3 He, 4 He, 
and 7 Li. 

2. 1 Theory 

Conditions for the synthesis of the light elements were attained in the early 
Universe at temperatures T > 1 MeV. In the early Universe, the energy 
density was dominated by radiation with 

9= — (2 + - + -AUT 4 (1) 

from the contributions of photons, electrons and positrons, and N v neutrino 
flavors (at higher temperatures, other particle degrees of freedom should be 
included as well). At these temperatures, weak interaction rates were in equi- 
librium. In particular, the processes 

n + e + <-> p + D e 
n + v e <-> p + e~ 

n <-* p + e~ + v e (2) 

fix the ratio of number densities of neutrons to protons. At T ^> 1 MeV, 
(n/p) ~ 1. 

The weak interactions do not remain in equilibrium at lower temper- 
atures. Freeze-out occurs when the weak interaction rate, T w k ~ G^T 5 
falls below the expansion rate which is given by the Hubble parameter, 
H - VUFp - T 2 /M P , where M P = 1/VGn ^ 1.2 x 10 19 GeV. The f3- 
interactions in eq. (||) freeze-out at about 0.8 MeV. As the temperature falls 
and approaches the point where the weak interaction rates are no longer fast 
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enough to maintain equilibrium, the neutron to proton ratio is given approx- 
imately by the Boltzmann factor, (n/p) ~ e - Am / T ^ 1/6, where Am is the 
neutron-proton mass difference. After freeze-out, free neutron decays drop 
the ratio slightly to about 1/7 before nucleosynthesis begins. 

The nucleosynthesis chain begins with the formation of deuterium by 
the process, p + n — > D + 7. However, because of the large number of 
photons relative to nucleons, rj^ 1 — n^/ns ~ 10 10 , deuterium production is 
delayed past the point where the temperature has fallen below the deuterium 
binding energy, Eb = 2.2 MeV (the average photon energy in a blackbody 
is E-y ~ 2. IT). This is because there are many photons in the exponential 
tail of the photon energy distribution with energies E > Eb despite the fact 
that the temperature or E 1 is less than Eb- The degree to which deuterium 
production is delayed can be found by comparing the qualitative expressions 
for the deuterium production and destruction rates, 

T p w n B crv (3) 
T d w n 7 ave~ EB/T 

When the quantity r)~ 1 exp(—EB /T) <~ 1, the rate for deuterium destruction 
(D + 7 — > p + n) finally falls below the deuterium production rate and the 
nuclear chain begins at a temperature T ~ O.lMeV. 

In addition to the p (n, 7) D reaction, the other major reactions leading 
to the production of the light elements are: 

D (D, p) T D (n, 7 ) T 3 Hc (n, p) T 

D (D, n) 3 He D(p, 7 ) 3 He 
Followed by the reactions producing 4 He: 
D (D, 7) 4 Hc 3 He ( 3 He, 2p) 4 He 

D ( 3 He, p) 4 He T(p, 7 ) 4 Hc 

T (D, n) 4 He 3 Hc(n, 7 ) 4 He 
The gap at A = 5 is overcome and the production of 7 Li proceeds through: 
3 He ( 4 He,7) 7 Be 

-» 7 Li + e+ + v e 

T ( 4 He, 7 ) 7 Li 
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Figure 1. The nuclear network used in BBN calculations. 



The gap at A = 8 prevents the production of other isotopes in .any significant 
quantity. The nuclear chain in BBN calculations was extendeco and is shown 
in Figure |l|. 

The dominant product of big bang nucleosynthesis is 4 He and its abun- 
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dance is very sensitive to the (n/p) ratio 

Yp = ■ 2 ^ - „ 0.25 (4) 
P [1 + (n/p)] V ' 

i.e., an abundance of close to 25% by mass. Lesser amounts of the other light 
elements are produced: D and 3 He at the level of about 10 -5 by number, and 
7 Li at the level of 10 -10 by number. 

Historically, BBN as a theory explaining the observed element abun- 
dances was nearly abandoned due its inability to explain all element abun- 
dances. Subsequently, stellar nucleosynthesis became the leading theory for 
element productionEl. However, two key questions persisted. 1) The abun- 
dance of 4 He as a function of metallicity is nearly flat and no abundances 
are observed to be below about 23% as exaggerated in Fig. ||. In partic- 
ular, even in systems in which an element such as Oxygen, which traces 
stellar activity, is observed at extremely low values (compared with the so- 
lar value of O/H as 8.5 x 1CP 4 ), the 4 He abundance is nearly constant. 
This is very different from all other element abundances (with the excep- 
tion of 7 Li as we will see below). For example, in Figure [|, the N/H 
vs. O/H correlation is shown. As one can clearly see, the abundance of 
N/H goes to 0, as O/H goes to 0, indicating a stellar source for Nitrogen. 
2) Stellar sources can not produce the observed abundance of D/H. Indeed, 
stars destroy deuterium and no astrophysical site is known for the production 
of significant amounts of deuteriumQ. Thus we are led back to BBN for the 
origins of D, 3 He, 4 He, and 7 Li. 



« M 4 He is Primordial! 
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Figure 2. The 4 He mass fraction as determined in extragalactic H II regions as a function 
of O/H. 
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Figure 3. The Nitrogen and Oxygen abundances in the same extragalactic HII regions with 
observed 4 He shown in Figure | 



The resulting abundances of the light elements are shown in Figure [|, 
over the range in 7710 = 10 10 ?7 between 1 and 10. The left plot shows the 
abundance of 4 He by mass, Y, and the abundances of the other three isotopes 
by number. The curves indicate the central predictions from BBN, while the 
bands correspond to the uncertainty in the predicted abundances based pri- 
marily the uncertainty in the input nuclear reactions as computed by Monte 
Carlo in ref. El. This theoretical uncertainty is shown explicitly in the right 
panel as a function of 7710- The dark shaded boxes correspond to the observed 
abundances of 4 He and 7 Li and will be discussed below. The dashed boxes 
correspond to the ranges of the elements consistent with the systematic un- 
certainties in the observations. The broad band shows a liberal range for 7710 
consistent with the observations. 

At present, there is a general concordance between the theoretical pre- 
dictions and the observational data, particularly, for 4 He and 7 Li. These two 
elements indicate that 77 lies in the range 1.7 < 77 < 4.7, corresponding to a 
range in fi^/i 2 = 0.006 — 0.0170. There is limited agreement for D/H as well, 
as will be discussed below. D/H is compatible with 4 He and 7 Li at the 2a 
level in the range 4.7 < 77 < 6.2 (fl B h 2 = 0.017 - 0.023). 



a f2 is the total density of matter relative to the critical density and f!g is the fraction of 
critical density in baryons. h is Hubble parameter scaled to 100 km Mpc — 1 s — 1 . 
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Figure 4. The light element abundances from big bang nucleosynthesis as a function of rjio. 



2.2 Data- 4 He 

The primordial 4 Hc abundance is best determined from observations of Hell — > 
Hel recombination lines in extragalactic HII (ionized hydrogen) regions. There 
is a good collection of abundance information on the 4 He mass fraction, Y, 
O/H, and N/H in over 70 such regionsa'EI. Since 4 He is produced in stars along 
with heavier elements such as Oxygen, it is then expected that the primordial 
abundance of 4 He can be determined from the intercept of the correlation 
between Y and O/H, namely Y p = Y(0/H — ► 0). A detailed analysialj of the 
data found 



Yn = 0.238 ± 0.002 ±0.005 



(5) 



The first uncertainty is purely statistical and the second uncertainty is an 
estimate of the systematic uncertainty in the primordial abundance determi- 
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Figure 5. The Helium (Y) and Oxygen (O/H) abundances in extragalactic HII regions, 
from refs. □ and □. Lines connect the same regions observed by different groups. The 
regression shown leads to the primordial 4 He abundance given in Eq. (51). 



nation. The solid box for 4 He in Figure [| represents the range (at 2cr sta t) 
from (|^). The dashed box extends this by including the systematic uncer- 
tainty. The 4 He data is shown in Figure la. Here one sees the correlation of 
4 He with O/H and the linear regression which leads to primordial abundance 
given in Eq. (|J). 

The helium abundance used to derive (||) was determined using assumed 
electron densities n in the HII regions obtained from SII data. Izotov, Thuan, 
& LipovetskyEl proposed a method based on several He emission lines to "self- 
consistently" determine the electron density. Their data using this method 
yields a higher primordial value 

Y p = 0.244 ± 0.002 ± 0.005 (6) 

One should also note that a recent determination!! 4 ! of the 4 He abundance 
in a single object (the SMC) also using the self consistent method gives a 
primordial abundance of 0.234 ± 0.003 (actually, they observe Y = 0.240 ± 
0.002 at [O/H] = -0.8, where the abundance [O/H] is defined as the log of 
the abundance relative to the solar abundance,JX/H] = log([X/H]/[X/H] Q )). 
Indeed, this work was extended in a reanalysisE] of selected regions from ref. 
□ and argue for corrections which lower the abundance to 0.239 ±0.002 based 
on a fit to 5 regions. 

As one can see, the resulting primordial 4 He abundance shows significant 
sensitivity to the method of abundance determination, leading one to con- 
clude that the systematic uncertainty (which is already dominant) may be 
underestimatedEj. 
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2.3 Data- 7 Li 

The abundance of 7 Li has been determined by observations of over 100 hot, 
population-II stars, and is found to have a very nearly uniform abundanceEj. 
For stars with a surface temperature T > 5500 K and a metallicity less than 
about l/20th solar (so that effects such as stellar convection may not be 
important), the abundances show little or no dispersion beyond that which is 
consistent with the errors of individual measurements. There is, however, an 
important source of systematic error due to the possibility that Li has been 
depleted in these stars, though the lack of dispersion in the Li data limits 
the amount of depletion. In fact, a small observed slope in Li vs Fe and the 
tiny dispersion about that correlation indicates that depletion is negligible in 
these stars Furthermore, the slope may indicate a lower abundance of Li 
than that in (6). The observation^ of the fragile isotopeJjLi is another good 
indication that 7 Li has not been destroyed in these staralH. _ 

The weighted mean of the 7 Li abundance in the sample of ref. tB is [Li] = 
2.12 ([Li] = log 7 Li/H + 12). It is common to test for the presence of a slope in 
the Li data by fitting a regression of the form [Li] = a + (3 [Fe/H]. These data 
indicate a rather large slope, (3 = 0.07 — 0.16 and hence a downward shift 

in the "primordial" lithium abundance A [Li] = —0.20 0.09. Models of 

galactic evolution which predict a smalLslope for [Li] vs. [Fe/H], can produce 
a value for (3 in the range 0.04 - 0.07 El Overall, when the regression based 
on the data and other systematic effects are taken into account a best value 
for Li/H was found to bala 

Li/H = 1.23 ±0.1 x 10" 10 (7) 

with a plausible range between 0.9 - 1.9 xl0~ 10 . The dashed box in Figure 
H corresponds to this range in Li/H. 

Figure H shows the different Li components for a model with ( 7 Li/H) p = 
1.23 x 10 -10 . The linear slope produced by the model is independent of the 
input primordial value (unlike the log slope given above) . The model of ref. lij 
includes in addition to primordial 7 Li, lithium produced in galactic cosmic ray 
nucleosynthesis (primarily a + a fusion), and Li produced by the ^-process 
during type II supernovae. As one can see, these processes are not sufficient 
to reproduce the population I abundance of 7 Li, and additional production 
sources are needed. 

2.4 Likelihood Analyses 

At this point, having established the primordial abundance of at least two 
of the light elements, 4 He and 7 Li, with reasonable certainty, it is possible 
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Figure 6. Confernbutions to the total predicted lithium abundancp-from the adopted GCE 
model of ref. £3, compared with low metallicity stars (from ref. 113) and a sample of high 
metallicity stars. The solid curve is the sum of all components. 



to test the concordance of BBN theory with observations. Two elements are 
sufficient for not only constraining the one parameter theory of BBN, but also 
for testing for consistency^]. A theoretical likelihood function for 4 He can be 
defined as 

L BBN (Y, F BBN ) = e -(Y-r^(v)f/^ (8) 

where !bbn(??) is the central value for the 4 He mass fraction produced in 
the big bang as predicted by the theory at a given value of 77. cviif the 
uncertainty in that value derived from the Monte Carlo calculations &Q and 
is a measure of the theoretical uncertainty in the BBN calculation. Similarly 
one can write down an expression for the observational likelihood function. 
Assuming Gaussian errors, the likelihood function for the observations would 
take a form similar to that in (||). 

A total likelihood function for each value of r\ is derived by convolving 
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the theoretical and observational distributions, which for 4 He is given by 



L Ho totai(ry) 



dYL 



BBN 



(I^bbn (r,)) L (Y,Y ) 



(9) 



An analogous calculation is performed^ for 7 Li. The resulting likelihood 
functions from the observed abundances given in Eqs. (||) and ([?]) is shown in 
Figure 0. As one can see there is very good agreement between 4 He and 7 Li 
in the range of 7710 ~ 1.5 - 5.0. 




Figure 7. Likelihood distribution for each of 4 He and 7 Li, shown as a function of rj in the 
upper panel. The lower panel shows the combined likelihood function. 



The combined likelihood, for fitting both elements simultaneously, is given 
by the product of the two functions in the upper panel of Figure ^ and is shown 
in the lower panel. The 95% CL region covers the range 1.7 < 7710 < 4.7, with 
the peak value occurring at 7710 = 2.4. This range corresponds to values of 
ilsh 2 between 

0.006 < fl B h 2 < 0.017 (10) 
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with a central value of fish 2 = 0.009. Using the higher value for 4 He (Y p = 
0.244) would result in an upward shift in 7710 by about 0.1 - 0.2. 

2.5 Data-D/H 

The remaining two light elements produced by BBN are D and 3 Hc. For 
the most part, the abundances of these elements are determined either in the 
local interstellar medium or in our own solar system. As such, one needs 
a model of galactic chemical evolution to tie them to the BBN abundances. 
Deuterium is predicted to be a monotonically decreasing function of time. 
The degree to which D is destroyed, however, is a model dependent question 
and related to the production of 3 Hc. Stellar models predict-that substantial 
amounts of 3 He are produced in stars between 1 and 3 M Q El. It should be 
emphasized that this predictionJa in fact consistent with the observation of 
high 3 He/H in planetary nebulaeE3. However, the implementation of standard 
model 3 He yields in chemical evolution_models leads. to an overproduction of 
3 Hc/H particularly at the solar epochEl See ref. Ej for attempts to resolve 
this problem. Because of this model dependence, I will not consider 3 He 
further here. 

Despite the problem of relating many of the locally observed D /H mea- 
surements to BBN, we can use the ISM values to set a firm lower limit on 
primordial D/H, and hence an upper limit to 77 because of the monotoni- 
cally decreasing history of D /H in the Galaxy. This is shown as the dashed 
(half)-box in Fig. |. 

There have been several reported measurements of D/H in high redshift 
quasar absorption systems. Such measurements are in principle capable of 
determining the primordial value for D/H and hence 77, because of the strong 
and monotonic dependence of D/H on 77. However, at present, detections 
of D/H using quasar absorption systems do not yield a conclusive value for 
D/H. In addition to the earlier determinationsE3 in two Lyman limit systems 
of 4.0i_65 and 3.25 ± 0.3 x 10~ 5 and the upper limit in a third Lyman limit 
systemt3 of 6.8 xl0~ 5 , there have been_three new measurements of D/HJjl 
Damped Lyman a systems: 2.5 ± 0.25 Ell; 2.25 ± 0.65 N; and 1.65 ± 0.35 E3; 
all times 10~ 5 . This data is shown in Fig. ||. Also shown in Fig. ||, are 
the adjusted D/H abundances when more complex velocity distributions are 
included as in refs. El 

Is there a real dispersion in D/H in these high redshift systems? There 
are two possible trends that go in a similar direction. The data may show an 
inverse correlation of D/H abundance with Si ErE£l. This may be an artifact of 
poorly determined Si abundances, or (as yet unknown) systematics affecting 
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Figure 8. The D/H data as a function of metallicity given by [Si/H] 



the D/H determination in high-column density (damped Lyman-a, hereafter 
DLA) or low-column density (Lyman limit systems) absorbers. On the other 
hand, if the correlation is real it would indicate that chemical evolution pro- 
cesses have occurred in these systems. The second trend is that the data may 
show an inverse correlation of D/H abundance with HI column density. If 
real, this would suggest that in the high column density DLA systems, which 
are most likely to have undergone some star formation, some processing of 
D/H must similarly have occurred at high redshift. One can only conclude: 
if the dispersion in D/H is real, it has profound consequences, as it indicates 
that some processing of D/H must, have occurred even at high redshift. 

It is interesting to speculated that the possible high redshift destruction 
of D/H is related to recent observationa_which suggest the existence of a white 
dwarf population in the Galactic halc£3. These observations could be signa- 
tures of an early population of intermediate-mass stars. Such a population 
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requires a Population III initial mass function different from that of the solar 
neighborhood. Also, to avoid overproduction of C and N, it is required that 
the Z = yields of these stars have low (~ 10~ 3 solar) abundances as sug- 
gested by some recent calculations, llnder these assumptions, it is possible 
to model the observed D vs Si trendEil. Such a scenario predicts a high cos- 
mic Type la supernova rate, while producing a white dwarf population that 
accounts for only ~ 1.5% of the dark halo. 

It is clear that a simple average of D /H abundance determinations does 
not make sense, at least without a proper enlargement of the error in the mean 
due to the poor \ 2 that such a mean would produce. Moreover, if deuterium 
destruction has occurred, we must question the extent to which any of these 
systems determine the value of Qb- It is important to note however, that for 
the upper end of the range (~ 5 x 10~ 5 ) shown in Fig. ||, all of the element 
abundances are consistent as will be discussed below. 

2.6 More Analysis 

It is interesting to compare the results from the likelihood functions of 4 He and 
7 Li with that of D/H. To include D/H, one would proceed in much the same 
way as with the other two light elements. We compute likelihood functions 
for the BBN predictions as in Eq. (||) and the likelihood function for the 
observations. These are then convolved as in Eq. (||). 

Using D/H = (3.0 ±0.3) x 10~ 5 as the primordial abundance, one obtains 
the likelihood function shown (shaded) in the upper panel of Fig. |[ The 95% 
CL region covers the range 5.0 < rjio < 7.4, with the peak value occurring at 
7/io = 5.8. This range corresponds to values of fls between 

0.018 < n B h 2 < 0.027 (11) 

with a central value of Fish 2 — 0.021. 

The combined likelihood, for fitting all elements (D, 4 He, and 7 Li) simul- 
taneously, is shown by the shaded curve in the lower panel of Figure Note 
that it has been scaled upwards by a factor of 25. In this case, the 95% CL 
region covers the range 4.7 < 7710 < 6.2, with the peak value occurring at 
7710 = 5.3. This range corresponds to values of Qb between 

0.017 < n B h 2 < 0.023 (12) 

with a central value of ilsh 2 = 0.019. 

It is important to recall however, that the true uncertainty in the low D/H 
systems might be somewhat larger. If we allow D/H to be as large as 5 x 10 -5 , 
the peak of the D/H likelihood function shifts down to rjio ~ 4. In this case, 
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Figure 9. Likelihood distribution for D/H, shown (shaded) as a function of r\ in the upper 
panel. Also shown are the 4 He and 7 Li likelihoods from Fig. Q The lower panel shows the 
combined likelihood function (shaded) compared to the previous case neglecting D/H. 



there would be a near perfect overlap with the high 77 7 Li peak and since the 
4 He distribution function is very broad, this would be a highly compatible 
solution. Given our discussion in the previous section concerning the current 
status of the D/H data, it is premature to claim a lack of concordance between 
BBN theory and observations. 



3 The CMB - BBN connection 

It is interesting to note the role of BBN in the prediction of the microwave 
background!! The argument is rather simple. BBN requires tempera- 
tures greater than 100 keV, which according to the standard model time- 
temperature relation, UT^y — 2.4/ V~N, where N is the number of relativistic 
degrees of freedom at temperature T, and corresponds to timescales less than 
about 200 s. The typical cross section for the first link in the nucleosynthetic 
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chain is 

av{p + n^ D + j) ~5x l(r 20 cm 3 /s (13) 
This implies that it was necessary to achieve a density 

n _ ^ 10 17 cm" 3 (14) 

avt v ; 

The density in baryons today is known approximately from the density of 
visible matter to be tib ~ 10~ 7 cm -3 and since we know that that the 
density n scales as i?~ 3 ~ T 3 , the temperature today must be 

T = (n Bo /n) 1/a TBBS ~ 10K (15) 
thus linking two of the most important tests of the Big Bang theory 




Figure 10. A\ 2 calculated with the MAXIMA-1 and COBE data as a function of parameter 
value. Solid blue circles show grid points in parameter space, and the green lines were 
obtained by interpolating between grid points. The parameter values where the green line 
intercepts the red dashed (dotted) line corresponds to the 68% (95%) frequentist confidence 
regionLl 
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Microwave background anisotropy measurements. Have made tremendous 
advances in the last few years. The power spectrurncJtiSEa has been measured 
relatively accurately out to multipolc moments corresponding to t ~ 1000. 



The details of this spectrum enable one to make accu 
large number of fundamental cosmological parametcrsLpll 
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of these results as found by a recent frequentist analysis£3 is shown in Fig. |lO| . 

The CMB anisotropies thus independently test the BBN prediction of 
Qsh 2 . At present, the predicted BBN baryon densities from D/H agree to 
an uncannitJevel with the most recent CMB resultsoEl The recent result 
from DASH2 indicates that fl B h 2 = 0.022±g:gg|, while that of BOOMERanG- 
98 0, fish 2 = 0.02lto'oo3 (using la errors) which should be compared to the 
BBN prediction given ixLeq. [ll]. These determinations are lower than value 
found by MAXIMA-1 B which yields VL B h 2 = 0.0325 ± 0.006. Given the 
current uncertainties, these results are consistent as can be seen in Fig. ITT] 
based on the recent frequentist analysis^ which found fish 2 = 0.026to' 06 
In addition, the BOOMERanG and DASI determinations are higher than the 
value fish 2 — 0.009 based oni-He and 7 Li. However, the measurements of 
the Cosmic Background ImageiH at smaller angular scales (higher multipoles) 
agree with lower BBN predictions and claims a maximum likelihood value for 
Sigh 2 = 0.009 (albeit with a large uncertainty). 
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Figure 11. Two-dimensional frequentist confidence regions in the (H , Sin) planeEJ The 
red, orange and yellow regions correspond to the 68%, 95%, and 99% confidence regions 
respectively. Standard calculations from big bang nucleosynthesis and observations of ^ 
predict a 95% confidence region of Sl B h 2 = 0.021 to 003 ' indicated by the shaded region. 
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4 Constraints from BBN 



Limits on particle physics beyond the standard model arc mostly sensitive to 
the bounds imposed on the 4 He abundance. As discussed earlier, the neutron- 
to-proton ratio is fixed by its equilibrium value at the freeze-out of the weak 
interaction rates at a temperature Tf ~ 1 MeV modulo the occasional free 
neutron decay. Furthermore, freeze-out is determined by the competition 
between the weak interaction rates and the expansion rate of the Universe 

G F 2 T/ ~ r wk (T/) = H(T f ) ~ ^G N NT f 2 (16) 

where N counts the total (equivalent) number of relativistic particle species. 
At T ~ 1 MeV, N = 43/4. The presence of additional neutrino flavors (or any 
other relativistic species) at the time of nucleosynthesis increases the overall 
energy density of the Universe and hence the expansion rate leading to a larger 
value of Tf, (n/p), and ultimately Yp. .-Because of the form of Eq. (|l6| ) it is 
clear that just as one can place limitscJ on N, any changes in the weak or 
gravitational coupling constants can be similarly constrained (for a discussion 
see ref. ell). 

Changes in N u actually affect not only 4 He, but also the abundances of the 
other light elements as seen in Fig. Hb3. We see the typical large dependence 
on N v in 4 He, but also note the shifts in the other elements, particularly D, 
and also Li over some ranges in r\. Because of these variations, one is not 
restricted to only 4 He in testing N v and particle physics. 

Just as 4 He and 7 Li wey .sufficient to determine a value for rj, a limit on 
N v can be obtained as wellE3c3o. The likelihood approach utilized above 
can be extended to include N v as a free parameter. Since the light element 
abundances can be computed as functions of both r\ and N„, the likelihood 
function can be defined bylfil replacing the quantity Ybbn {v) m ec l- @ with 
Ybbn (>7, N„) to obtain L He totai(?7, N v ). Again, similar expressions are needed 
for 7 Li and D. _ 

The likelihood distribution derived from the analysis of ref. c3 is shown 
in Fig. [l^ where iso-likclihood contours representing 68, 95, and 99 % CL 
are projected onto the 7710 — N v plane. As one can see, when only 4 He and 
7 Li are used, the allowed range in 7710 is rather broad and the upper limit 
to N„ <> 4.2 at relatively low 77. At present, one can not use D/H to fix N v , 
unless one has additional information on the value of 77, e.g., from the CMB as 
is demonstrated in Fig. [l4|. On the other hand when all three light elements 
are used, the narrow range in rj at relatively high values, allows for a tight 
constraint on N v < at 95 % CL. Furthermore, it is clear that the most realistic 
value, N v = 3 is in fact consistent with all of the observed abundances at the 
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Figure 12. BBN abundance predictionsEJ as a function of the baryon-to-photon ratio rj, for 
Nu = 2 to 7. The bands show the la error bars. Note that for the isotopes other than Li, 
the error bands are comparable in width to the thickness of the abundance curve shown. 
All bands are centered on N v = 3. 



2 a level. 

Of course as discussed above, the CMB provides an independent deter- 
mination of r], we can use that to fix rj, and test the available constraints on 
N y l3. For example, we can fix r\ = 5.8 x 1CP 10 which corresponds to the 
BOOMERanG and DASI determinations, and use the current 4 He and 7 Li 
data to compute a BBN likelihood function for N v . This is shown in the first 
panel of Fig. |TJ. Assuming a 10% uncertainty in 77, the 95 % CL range on 



Otai: submitted to World Scientific on February 1, 2008 19 



95% CL 
99% CL 




Figure 13. (a) Likelihood contours representing 68, 95, and 99 % CL projected onto the 
VlO — N v plane. In a) only 4 He and 7 Li are used. In b) only D is used, and in c) all 
three light elements are used. The crosses denote the position of the peak of the likelihood 
distribution. 



N v is 1.8 - 3.3, whereas the 1 sided upper limit with a N v > 3 priorCJ is 3.5. 
Using instead the D/H data, at the same baryon density, one finds an upper 
limit N v < 6.3. This likelihood is shown in panel b) of Fig. O. Alternatively, 
we could choose a lower value of rj — 2.4 x 10~ 5 , and using^He and 7 Li we 
find that the 95 % CL range is 2.2 - 3.9. This case is exemplified in panel c). 

5 Something's Missing 

There is considerable evidence for dark matter in the UniverseS. The best 
observational evidence is found on the scale of galactic halos and comes from 
the observed flat rotation curves of galaxies. There is also good evidence for 
dark matter in elliptical galaxies, as well as clusters of galaxies coming from X- 
ray observations of these objects. In theory, we expect dark matter because 1) 
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Figure 14. (a) The distribution irLiV„ assuming a value of r\ = 5.8 X 10 and primordial 
4 He & 7 Li abundances as in eqs. H and M. (b) As in (a) assuming a D measurement at 
the current precision (taken to be 3.0 ± 0.3 X 10 ). (c) As in (a), but with a value for 
r? = 2.4 X 10" 10 . 



inflation predicts Q = 1, and the upper limit on the baryon (visible) density of 
the Universe from big bang nucleosynthesis is Hb < 0.1; 2) Even in the absence 
of inflation (which does not distinguish between matter and a cosmological 
constant), the large scale power spectrum is consistent with a cosmological 
matter density of O ~ 0.3, still far above the limit from nucleosynthesis; 
and 3) our current understanding of galaxy formation is inconsistent with 
observations if the Universe is dominated by baryons. 

Indeed, we now have direct evidence from CMB anisotropy measurements 
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Figure 15. Two-dimensional frequentist confidence regions in the (Hm^a) plane. The 
red, orange and yellow regions correspond to the 68%, 95%, and 99% confidence regions 
respectively. The dashed black line corresponds to a flat universe, Q = Om + £2 A = !• 



that fitot — 1 to with in about 20 %. The CMB anisotropy allows for a 
determination of the curvature of the Universe. Therefore, while one can 
not determine unambiguously the value of fi ma tter; one can nx the sum of 
the matter and dark energy (cosmological constant) contributions. This is 
shown in Fig. [l5| from the frequentist analysis ir£3. When combined with high 
redshift sunernovae measurements indicating that the Universe is currently 
acceleratingcZl, one would conclude that the dark energy constitutes about 
65 % of the closure density leaving the remaining 35 % for matter. Given 
that the baryon density is less than 5 %, this leaves us with about 30 % of 
closure density for non-baryonic dark matter. Other dynamical arguments 
estimating the density of matter are in agreement with this allotment (see 
ref. c3 for a recent review). 

In fact, there are many reasons why most of the dark matter must be non- 
baryonic. In addition to the problems with baryonic dark matter associated 
with nucleosynthesis or the growth of density perturbations, it is very difficult 
to hide baryons. There are indeed very good constraints on the possible forms 
for baryonic dark matter in our galaxy. Strong cases can be made against hot 
gas, dust, jupiter size objects, and stellar remnants such as white dwarfs and 



In what follows, I will focus on supersymmetric candidates in which the 
relic abundance of dark matter contributes a significant though not excessive 
amount to the overall energy density. Denoting by Cl x the fraction of the 
critical energy density provided by the dark matter, the density of interest 
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falls in the range 

0.1 < n x h 2 < 0.3 (17) 

The lower limit in eq.(^) is motivated by astrophysical relevance. For lower 
values of Q x h 2 , there is not enough dark matter to play a significant role in 
structure formation, or constitute a large fraction of the critical density. The 



upper bound in (17), on the other hand, is an absolute constraint, derivable 



from the age of the Universe, which can be expressed as 

H t = [ dx(i-n-n A + n A x 2 + n/ x y 1/2 (is) 

Jo 

Given a lower bound on the age of the Universe, one can establish an upper 
bound on Qh 2 from eq.(|l^). The limit to > 12 Gyr translates into the upper 
bound given in (|l7|). Adding a cosmological constant does not relax the 
upper bound on (Tn 2 , so long as O + £1 a < 1 . If indeed, the_Ladications for a 
cosmological constant from recent supernovae observationscH turn out to be 
correct, the density of dark matter will be constrained to the lower end of the 
range in (p~7|) . Indeed for h 2 ~ 1/2, we expect il x h 2 ~ 0.15. 

6 Supersymmetric Dark Matter 

Although there are many reasons for considering supersymmetry as a can- 
didate extension to the standard model of strong, weak and electromagnetic 
interactionsE3, ope of the most compelling is its role in understanding the hi- 
erarchy problemtJ namely, why/how is m\y <C Mp. One might think naively 
that it would be sufficient to set mw <C Mp by hand. However, radiative 
corrections tend to destroy this hierarchy. For example, one-loop diagrams 
generate 



Sm 2 w = 



A 2 » m 2 w (19) 



where A is a cut-off representing the appearance of new physics, and the in- 
equality in ( |l9| ) applies if A ~ 10 3 TeV, and even more so if A ~ vtlqut ~ 10 16 
GeV or ~ Mp ~ 10 19 GeV. If the radiative corrections to a physical quan- 
tity are much larger than its measured values, obtaining the latter requires 
strong cancellations, which in general require fine tuning of the bare input 
parameters. However, the necessary cancellations are natural in supersym- 
metry, where one has equal numbers of bosons B and fermions F with equal 
couplings, so that ( |l9| ) is replaced by 

8m 2 w = o(?-) \m 2 B -m 2 F \ . (20) 
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The residual radiative correction is naturally small if \m% — m 2 F \ < 1 TeV 2 . 

In order to justify the absence of superpotential terms which can be re- 
sponsible for extremely rapid proton decay it is common in the minimal super- 
symmetric standard model (MSSM) to assume the conservation of R-parity. 
If R-parity, which distinguishes between "normal" matter and the supersym- 
metric partners and can be defined in terms of baryon, lepton and spin as 
R = (— i) 3B + L + 2S j i s unbroken, there is at least one supersymmetric particle 
(the lightest supersymmetric particle or LSP) which must be stable. 

The stability ,xif the LSP almost certainly renders it a neutral weakly in- 
teracting particleEx Strong and electromagnetically interacting LSPs would 
become bound with normal matter forming anomalously heavy isotopes. In- 
deed, there are very strong_Lipper limits on the abundances, relative to hy- 
drogen, of nuclear isotopesE^I, n/njj ^ 10~ 15 to 10~ 29 for 1 GeV <^ m <, 
1 TeV. A strongly interacting stable relics is expected to have an abundance 
n/n,H 10~ 10 with a higher abundance for charged particles. 

There are relatively few supersymmetric candidates which are not col- 
ored and are electrically neutral. The sneutrinco is one possibility, hiit in 
the MSSM, it has been excluded as a dark matter candidate by directEJ and 
indirectO searches. In fact, one can set an accelerator based limit on the sneu- 
trino mass from neutrino counting, mp > 43 GeV Ell. In this case, the direct 
relic .searches in underground low-background experiments require rriy > 1 
TeV Ea Another possibility is the gravitino which is probably the most diffi- 
cult to exclude. I will concentrate on the remaining possibility in the MSSM, 
namely the neutralinos. 



6.1 Relic Densities 

There are four neutralinos, each of which is a linear combination of the R = — 1 
neutral fermions,E2l: the wino W 3 , the partner of the 3rd component of the 
SU(2)l gauge boson; the bino, B, the partner of the U(l)y gauge boson; and 
the two neutral Higgsinos, Hi and H 2 . Assuming gaugino mass universality 
at the GUT scale, the identity and mass of the LSP are determined by the 
gaugino mass m 1 / 2 (or equivalently by the SU(2) gaugino mass M 2 at the 
weak scale - at the GUT scale, M 2 = m.1/2), the Higgs mixing mass fx, and 
the ratio of Higgs vevs, tan/3. In general, neutralinos can be expressed as a 
linear combination 



X = aB + f3W 3 + 7# x + SH 2 (21) 
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The solution for the coefficients a, f3, 7 and 5 for neutralinos that make up the 
LSP can be found by diagonalizing the mass matrix 



M 2 




V2 

1 32 V2 
\ V2 





Mi 

givi 
-giv2 

V2 



-92V1 

givi 




92l'2 

-giV2 
V2 

-A* 





B 
Hi 



(22) 



where Mi is a soft supersymmetry breaking term giving mass to the U(l) 
gaugino. In a unified theory Mi = M2 at the unification scale (at the weak 
scale, Mi = l^-Ma). As one can see, the coefficients a, (3, 7, and £ depend 
only on mi/2, Mi pjd tan/3. 

In Figure [l6]E3, regions in the M2, /Lt plane with tan/3 = 2 are shown in 
which the LSP is one of several nearly pure states, the photino, 7, the bino, 
B, a symmetric combination of the Higgsinos, -ff(i2) 5 or the Higgsino, S = 
sm(3Hi +cos/3i?2- The dashed lines show the LSP mass contours. The cross 
hatched regions correspond to parameters giving a charging (M /± ,i? ± ) state 
with mass < 45GeV ancLas such are excluded by LEPo. This constraint 
has been extended by LEPEd and is shown by the light shaded region and 
corresponds to regions where the chargino mass is < 104 GeV. The newer 
limit does not extend deep into the Higgsino region because of the degeneracy 
between the chargino and neutralino. Notice that the parameter space is 
dominated by the B or -H12 pure states and that the photino only occupies 
a small fraction of the parameter space, as does the Higgsino combination S. 
Both of these light states are now experimentally excluded. 

The relic abundance of LSP's is determined by solving the Boltzmann 
equation -for the LSP number density in an expanding Universe. The 
techniqueEa used is similar to that for computing the relic abundance of mas- 
sive neutrinos^. The relic density depends on additional parameters in the 
MSSM beyond M2,//, and tan/3. These include the sfermion masses, mj, 
the Higgs pseudo-scalar mass, m^, and the tri-linear masses A as well as two 
phases 9^ and 9a- To determine, the relic density it is necessary to obtain 
the general annihilation cross-section for neutralinos. In much of the param- 
eter space of interest, the LSP is a bino and the annihilation proceeds mainly 
through sfermion exchange as shown in Figure [T^. For binos, it is possible to 
adjust the sfermion masses to obtain closjare density in a wide mass-tange. Ad- 
justing the sfermion mixing parametersEJ or CP violating phasesE3l23 allows 
even greater freedom. 

Because of the p-wave suppression associated with Majorana fermions, 
the s-wave part of the annihilation cross-section is suppressed by the outgoing 
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Figure 16. Mass contours and composition of nearly pure LSP states in the MSSM 




Figure 17. Typical annihilation diagram for neutralinos through sfermion exchange. 



fermion masses. This means that it is necessary to expand the cross-section 
to include p-wave corrections which can be expressed as a term proportional 
to the temperature if neutralinos are in equilibrium. Unless the neutralino 
mass happens to lie near near a pole, such as m x ~ or m/j/2, in which 

case there are large contributions to the annihilation through direct s-channel 
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resonance exchange, the dominant contribution to the BB annihilation cross 
section comes from crossed i-channel sfermion exchange. In the absence of 
such a resonance, the thermally-averaged cross section for BB — + // takes 
the generic form 



v m\ ' 128tt 



B 



9 

m. 



+ {Yt + Y^){—^){l + ...)x 

= a + bx (23) 

where Yl(r) is the hypercharge of Jl(r), &f = ttA + — rnj, and we have 
shown only the leading P-wave contribution proportional to x = T/ms. Such 
an expansion yields very accurate results unless the LSP is near a threshold, 
near a resonance, or is nearly degenerate in mass with another SUSY particle 
in which cases a more accurate treatment is necessaryEll. Annihilations in 
the early Universe continue until the annihilation rate T ~ avn x drops below 
the expansion rate, H. For particles which annihilate through approximate 
weak scale interactions, this occurs when T ~ m x /20. Subsequently, the relic 
density of neutralinos is fixed relative to the number of relativistic particles. 

As noted above, the number density of neutralinos is tracked by a 
Boltzmann-likc equation, 

f = S^n-{av)(n 2 -nl) (24) 

where no is the equilibrium number density of neutralinos. By defining the 
quantity / = rt/T 3 , we can rewrite this equation in terms of x, as 



d f ( 1 2 2 



The solution to this equation at late times (small x) yields a constant value 
of /, so that n (x T 3 . The final relic density expressed as a fraction of the 
critical energy density can be written asc3 

where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon tem- 
perature with respect to x 5 due to the annihilations of particles with mass 
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to < Xfm x E3. The subscript / refers to values at freeze-out, i.e., when 
annihilations cease-. 

In Figure [l^ E2l, regions in the M2 — [i plane (rotated with respect to 
Figure [l6|) with tan/3 = 2, and with a relic abundance 0.1 < Qh 2 < 0.3 
are shaded. In Figure |l8|, the sfermion masses have been fixed such that 
too = 100 GeV (the dashed curves border the region when m = 1000 GeV). 
Clearly the MSSM offers sufficient room to solve the dark matter problem. In 
the higgsino sector P12, additional types of annihilation processes known as 
coannihilations EJOO between Hn 2 \ and the next lightest neutralino (H\i2]) 
must be included. These tend to significantly lower the relic abundance in 
much of this_sector and as one can see there is little room left for Higgsino 
dark matteiO. 




Figure 18. Regions in the M2— fi plane where 0.1 < flh 2 < 0.3 Ej- Also shown are the 
Higgsino purity contours (labeled 0.1, 0.5, and 0.9). As one can see, the shaded region is 
mostly gaugino (low Higgsino purity). Masses are in GeV. 
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As should be clear from Figures [l6| and binos are a good and likely 
choice for dark matter in the MSSM. For fixed m,-, Vlti 2 > 0.1, for all rng = 
20 — 250 GeV largely independent of tan /? and the sign of /i. In addition, the 
requirement that mr > nig translates into an upper bound of about 250 GeV 
on the bino mass@@. By further adjusting the trilinear A and accounting 
for sfermion mixing this upper bound can be relaxecO and by allowing for 
non-zero phases in the MSSM, the upper limit can be extended to about 600 
GeV lB. For fixed Qh 2 — 1/4, we would require sfermion masses of order 120 
- 250 GeV for binos with masses in the range 20 - 250 GeV. The Higgsino 
relic density, on the other hand, is largely independent of m^. For large 
fi, annihilations into W and Z pairs dominate, while for lower /i, it is the 
annihilations via Higgs scalars which dominate. Aside from a narrow region 
with m^ i2 < mw and very massive Higgsinos with rrijj 12 ft 500 GeV, the 
relic density of H12 is very low. Above about 1 TeV, these Higgsinos are also 
excluded. 

As discussed in section 4, one can make a further reduction in the number 
of parameters by setting all of the soft scalar masses equal at the GUT scale 
(similarly for the A parameters as well). In this case the free parameters are 

m o,nT-i/2,A and tan/3 , (27) 

with fi and being determined by the clcctroweak vacuum conditions, the 
former up to a sign. We_iefer to this scenario as the constrained MSSM 
(CMSSM). In Figure 19 £3, this parameter space is shown for tan/3 = 10. 
The light shaded region corresponds to the portion of parameter space where 
the relic density £l x h 2 is between 0.1 and 0.3. The dark shaded region (in the 
lower right) corresponds to the parameters where the LSP is not a neutralino 
but rather a f^. The cosmologically interesting region at the left of the 
figure is due to the appearance of pole effects. There, the LSP can annihilate 
through s-channel Z and h (the light Higgs) exchange, thereby allowing a very 
large value of mo. However, this region is excluded by the phenomenological 
constraints described below. The cosmological region extends toward large 
values of ™w«, due to coannihilation effects between the lightest neutralino 
and the tr L3. For._Don-zero values of A , there are new regions for which 
X — t are importantLJ. 

Figure |l9| also shows the current experimental constraints on the CMSSM 
parameter space provided by direct searches at LEP. One of-these is the limit 
m x ± > 103.5 GeV provided by chargino searches at LEP El LEP has also 
provided lower limits on slepton masses, of which the strongest is mg ft 99 
GeV 13. This is shown by dot-dashed curve in the lower left corner of Fig. |l^. 

Another important constraint is provided by the LEP lower limit on the 
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Figure 19. Compilation of phenomenological constraints on the CMSSM for tan (3 = 10, fi > 
0, assuming Ao = 0, m t = 175 GeV and m (j (m i) )|^ = 4.25 GeY L3. The near-vertical 
lines are the LERJimits m x ± = 103.5 GeV (dashed and black)Eil, and m h = 114.1 Ge¥ 
(dotted and red)Lfl. Also, in the lower left corner we show th mg = 99 GeV contounlj. 
In the dark (brick red) shaded regions, the LSP is the charged fi, so this region is ex- 
cluded. The ligkHjturquoisc) shaded areas are the cosmologically preferred regions wj%li 
0.1 < Qh 2 < 0.3L3. The medium (dark green) shaded regions are excluded by b — > s~/ 
The shaded (pink) regions in the upper right regions delineate the 2 <r range of — 2. The 
±1 cr contours are also shown as solid black lines inside the shaded region. 



Higgs mass: mn > 114.1 GeV This holds in the Standard Model, for 
the lightest Higgs boson h in the general MSSM for tan/3 < 8, and almost 
always in the CMSSM for all tan fi. Since rrih is sensitive to sparticle masses, 
particularly m^, via loop corrections, the Higgs limit also imposes important 
constraints on the CMSSM parameters, principally m^ji as seen in Fig. [H| 
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The constraints are evaluated using FeynHiggsEJ, which is estimated to have 
a residual uncertainty of a couple of GeV in m^. 

—Also shown in Fig. [l^ is the constraint imposed by measurements of b — ► 
sj 113. These agree with the Standard Model, and therefore provide bounds on 
MSSM particles, such as the chargino and charged Higgs masses, in particular. 
Typically, the b — > sj constraint is more important for (i < 0, but it is also 
relevant for [i > 0, particularly when tan/3 is large. The region excluded by 
the b — > sj constraint is the dark shaded (green) region to the left of the plot. 

The final experimental constraint considered is that due to the measure- 
ment of the anomalous magnetic moment of the muon. The BNL E821 exper- 
iment reported last year a new measurement of a M = \{g^ — 2) which deviated 
by 2.6 standard deviations from the best Standard Model prediction available 
at that time l3. The largest contribution to the errors in the comparison with 
theory was thought to be the statistical error of the experiment, which will 
soon be significantly reduced, as many more data have already been recorded. 
However, it has recently been realized that the sign of the most important 
pseudoscalar-meson pole part of the light-by-light scattering contributiorO 
to the Standard Model prediction should be reversed, which reduces the ap- 
parent experimental discrepancy to about 1.6 standard deviations. 

As many authors have pointed outEj, a discrepancy between theory and 
the BNL experiment could well be explained by supersymmetry. As seen in 
Fig. [l9[ this is particularly easy if fi > 0. The medium (pink) shaded region 
in the figure is the new 2 a allowed region: —6 < Sa^ x 10 10 < 58. With the 
change in sign of the meson-pole contributions to light-by-light scattering, 
good consistency is also possible for fi < so long as either mi/a or mo are 
taken sufficiently large. 

In addition to the coannihilation region discussed above, another mech- 
anism for extending the allowed CMSSM region to lajjgf m x is rapid anni- 
hilation via a direct-channel pole when m x ~ \vtla r tt 6 !. This may yield—a 
'funnel' extending to large m 1 / 2 and m at large tan/3, as seen in Fig. p0|O. 

6.2 Detection 

As an aid to the assessment of the prospects for detecting sparticles at dif- 
ferent accelerators, benchmark sets of supersymmetric parameters have often 
been found useful, since they provide a focus for_concentrated discussion. A 
set of proposed post-LEP benchmark scenariosEJ in the CMSSM are illus- 
trated schematically in Fig. [2l]. They take into account the direct searches for 
sparticles and Higgs bosons, b — > sj and the preferred cosmological density 
range. All but one of the benchmark points are consistent with — 2 at the 
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Figure 20. As in Fig. |l| for tan /3 = 50 EE 



2 a level. 

The proposed points were chosen not to provide an 'unbiased' statistical 
sampling of the CMSSM parameter space, but rather are intended to illustrate 
the different possibilities that are still allowed by the present constraintsE3. 
Five of the chosen points are in the 'bulk' region at small m 1 / 2 and mo, four 
are spread along the coannihilation 'tail' at larger mi/2 for various values of 
tan (3, two are in rapid-annihilation 'funnels' at large m.1/2 and mp. At large 



values of mo, (larger than that shown in Figs. [19 and 20, there is another 



region_where the cosmological range is satisfied, namely in the 'focus-point' 
regionE3 along the boundary where electroweak symmetry no longer occurs 
(shown in Fig. |2l] as the shaded region in the upper left corner.) Two points 
were chosen in the 'focus-point' region at large mo- The proposed points range 
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Figure 21. Schematic overview of the CMSSM benchmark points proposed in El The points 
are intended to illustrate the range of available possibilities. The labels correspond to the 
approximate positions of the benchmark points in the (mi/2 > mo) plane. They also span 
values of tan j3 from 5 to 50 and include points with fi < 0. 



. detection 
s5 ' for these 



over the allowed values of tan/? between 5 and 50. Prospects for t 
of supersymmetry at future colliders was studied extensively in 
points in particular. 

Because, the LSP as dark matter is present locally, there are many avenues 
for pursuing dark matter detection. Here I conclude by showing the prospects 
for direct detection for the benchmark points discussed aboveEH. Fig. ^ shopas 
rates for the elastic spin-independent scattering of_supersymmetric_relicsEll, 
including thejjrojected sensitivities for CDMS II E3 and CRESSIEJ (solid) 
and GENIUSEl (dashed). Also shown are the cross sections calculated in the 
proposed benchmark scenarios discussed in the previous section, which are 
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Figure 22. Elastic spin-independent scattering-qf supersymmetric relics on (a) protons and 
(b) neutrons qaiculated in benchmark scenariosEJ, compared with the projected sensitivities 
for CDMS II E3 and CRESSTE3 (soprf) and GENIUSE3 (dashed). The predictions of our 
code (blue crosses) and Neutdriveil23 (red circles) for neutralino-nucleon scattering are 
compared. The labels A, B, ...,L correspond to the benchmark points as shown in Fig. BU 



considerably below the DAMAo range (10 -5 — 1CP 6 pb), but may be within 
reach of future projects. Indirect searches for supersymmetric dark matter 
via the products of annihilations in the galactic halo or inside the Sun also 
have prospects in some of the benchmark scenariosEZl 
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